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INTRODUCTION 


This  report  summarizes  the  work  conducted  under  Contract  No. 

H0220025  during  the  period  1  March  -  1  September  1972.  The  program  Is  generally 
concerned  with  an  understanding  of  how  stress  waves  Interact  with  discrete 
faults,  randomly  distributed  faults  and  tunnels  In  rock  media.  The  study 
comprises  both  theoretical  analyses  and  two-d<menslonal  computer  calculations 
of  selected  Interaction  problems.  The  work  carried  out  during  this  period  was 
primarily  concerned  with  the  characterization  of  randomly  faulted  media  and 
the  development  of  an  algorithm  for  controlling  crack  propagation  In  the  SHEP 
code. 

2.  TECHNICAL  PROGRESS;  RANDOMLY  FAULTED  MEDIA 

2.1  POROSITY 

In  trying  to  understand  the  effects  of  strong  stress  waves  on  rock 

media,  e.g. ,  hard  granite.  It  Is  first  necessary  to  assign  a  length  scale 

within  which  damage  occurs.  This  length  scale  depends,  of  course,  on  the 

energy  associated  with  the  excavatloo  process  or  damage  process  caused  by 

a  high  explosive  Curst.  There  are  two  categories  of  bursts,  those  which  are 

used  for  excavation  (a  few  tons  of  high  explosive),  and  those  which  are 

designed  to  damage  underground  structures  or  tunnels  (several  kllotons  or 

more).  It  Is  generally  agreed^  that  the  scaled  elastic  radius  (the  radius  of 

a  spherical  region  within  which  plastic  yield  and  significant  fracture  damage 

1/3 

may  occur)  surrounding  a  burst  center  Is  of  the  order  of  400  ft/kt  '  .  This 
means  that  the  elastic  radius  associated  with  a  megaton  burst  Is  of  the  order 
of  4000  ft.  On  the  other  hand  the  elastic  radius  associated  with  an  excava¬ 
tion  burst  may  be  as  little  as  40  ft. 

From  tests  such  as  PILEDRIVER  and  HARD  HAT  which  were  conducted 

2 

at  Nevada  Test  Site  15  In  hard  granite.  It  Is  known  that  fractures,  faults  or 
blocks  are  distributed  with  spaclngs  significantly  less  than  1  ft.  Thus  a 
sphere  of  radius,  as  li;;tle  as  40  ft.  will  contain  a  large  number  of  flaws.  In 


1 


such  a  sphere  or  a  larger  one,  om;  Is  led  to  the  view  that  the  fracture  area 
is  so  complicated  and  so  diversified  that  one  may  profitably  approach  the 
problem  of  defining  the  flaw  structure  on  a  stochastic  basis,  and  assume  that 
the  flaw  area  Is  randomly  distributed.  This  Idea  has  been  subsumed  by  other 
Investigators^. 

With  this  thought  In  mind,  we  undertake  to  characterize  a  randomly 
faulted  medium  by  elastic  moduli  which  depend  on  the  fracture  area  per  unit 
vof.ume.  There  arise  Immediately  a  number  of  model  problems  that  one  can  use 
as  a  guide  to  this  characterization.  The  simplest  one  and  the  first  of  several 
we  shall  discuss  Is  the  hydrostatic  tension  of  a  porous  medium.  If  one 
makes  the  first  order  approximation  that  a  unit  cell  of  a  randomly  porous 
medium  Is  a  spherical  cell,  then  It  Is  easy  to  show  that  the  effective  bulk 
modulus  of  such  a  cell  Is  given  by 

1  -  V 

*^eff  ~  3^  (1) 

K  ^4G 


where 


Is  the  volume  fraction  of  pores,  and 
are  the  elastic  moduli  of  the  continuum 


It  Is  also  easy  to  relate  Vp  to  the  surface  area  per  unit  volume  o?,  namely 


''p=<t4 
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(2) 


where  b  Is  a  length  scale  which  has  to  do  with  porosity  or  flaw  size.  We 

suspect  that  interaction  among  the  adjacent  flaw  will  result  In  a  not  dissimilar 
relation  where: 

'‘p  =  (— )  (3) 


and  where  n  and  b  must  be  determined  experimentally. 
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Expression  (1)  is  based  on  linear  elastic  theory.  Because  the 
applied  stress  becomes  concentrated  at  the  pore  surface,  it  is  of  interest 
to  see  under  what  conditions(s)  (1)  breaks  down.  We  approach  the  question 
by  analyzing  the  effect  of  geometric  nonlinearity. 

Assume  that  a  spherical  shell  of  outer  radius,  b,  inner  radius,  a. 
Is  subjected  to  isotropic  tension,  S,  at  the  outer  surface,  i.e. 


a  1  =  +  S  (4) 

^  r=b 

We  desire  to  determine  the  amount  of  dilatation.  The  bulk  material  will 
actually  dilate  very  slightly,  while  the  pore  will  dilate  more  or  less  highly 
depending  upon  the  tension.  Let  us  assume  that  the  material  properties  are 
given  by; 

W  =  f  (LXj  -  3)  +Gfa)  (5) 

where  W  is  the  strain  energy  function 

is  the  stretch  ratio  in  the  1^^  principal  direction 
J  is  the  volume  ratio,  J  =  ffX^  (6) 

and  G  is  the  shear  modulus  of  the  material. 

The  Piola  stress  is  given  by: 


(11) 


Thus 


0  r 

dr 

dr 


We  can  also  write 


and 

Thus 


4  '• 


°e  .  .im 


(12) 

(13) 

(14) 

(15) 

(16) 


In  writing  ( 5  ) ,  we  have  included  in  the  first  term  on  the  right  hand 
side  all  the  geometrical  nonlinearity  associated  with  shear.  The  second  term 
has  been  left  arbitrary,  but  does  represent  the  effect  of  dilatation.  Since  J 
will  nowhere  throughout  the  spherical  shell  be  significantly  different  from 
unity,  it  matters  little  what  function  f(J)  we  choose  so  long  as  it  reduces  to 
Hooke's  law  in  the  limit  of  infinitesimal  strain.  When  we  make  this  identi¬ 
fication  ,  we  find  that; 

f(l)  =  0  (17) 

f'(l)=-l  (18) 


f"(l) 


1 _ K  +  I 

l-2v  ■  G  3 


G 


(19) 


Apart  from  these  mild  restrictions,  f  may  have  any  form  whatsoever.  We  are 
going  to  choose  f  to  make  the  ensuing  analysis  convenient. 

Elastostatic  equilibrium  is  guaranteed  by  the  statement: 


2  dr 


0  - 
r 


e 


=  0 


(20) 
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1 


We  also  have 


so  that: 


r  =  Xr 


and 


dr  -  ^  + 


r  -  l~ —  -  XI 
dr  “  '  ^  A\ 


dX 


(21) 


(22) 


(23) 


Equation  (22)  will  be  used  to  determine  r  as  a  function  of  X^  or 
Vice  versa,  when  J  Is  known, as  a  function  of  X.  Using  (23),  we  can  rewrite 
(20)  In  the  form: 


J 

<;;2 dT^  “ 


(24) 


After  Introducing  (15)  and  (16)  and  simplifying,  one  obtains 


^  .-2JXizJ),- 

X(l+X^f") 


(25) 


Equation  (24)  Is  Immediately  separable  If  we  set 


f"  = 


(26) 


Now  let 


so  that 


J  =  w^X^ 


dtnX  = - (C_+  w  )  dw 


w  [C  +1  w  +  jw^] 


(27) 


(28) 


Now  observe: 


asr-*b,X-l  +  c,J-'l+0(€^),w-'l-€ 
asr-a,  X-^,J  =  1  +  O(e^),  w -c 


(29) 

(30) 
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Thus,  at  either  limit,  the  governing  equation  behaves  like: 


d-tnX  »  -  d-t-nw 

(31) 

Thus,  to  an  excellent  approximation: 

O 

II 

(32) 

CO  «-• 

O 

II 

(0 

(3  3) 

Insertion  of  this  result  into  (23)  yields: 

3  ^2 

x=(c5  +  -t-) 

r 


The 


radial  component  of  the  Cauchy  stress  is  given  by: 


(34) 


^  x'"  J 

which,  with  (33)  ::nd  (34)  yields: 


^r  P  3  /Q  3  — 2_\ 
_  =  Cl  (Cj  +  3  ) 


-4/3 


1-  3C  - 


1_3C 
'  ^1 


(35) 


(36) 


We  now  set: 


0  =  S ,  at  r  =  b 

r 


0  =  0,  at  r  =  a 
r 


(37) 

(38) 


We  can  solve  (37)  and  (38)  to  an  excellent  approximation  by  first  setting 


=  1  +  0 


(39) 


and  then  linearizing  both  expressions  In  8.  This  T^XTole 

the  materials  o£  interest  have  Poisson's  ratio  v=  1/4,  or  C  =  2.  In 
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of  a  rubberlike  material,  this  approximation  breaks  down  because  C  is  more 
like  1000. 


If  we  define 


Xj  -  1  -  6 

then  we  arrive  at: 


B  = 


(40) 


(41) 


which  is  identical  with  (1) .  Thus  geometrical  nonlinearity  has  no  overriding 
effect  even  though  strain  at  the  pore  may  be  very  large  and  even  though  the 
pore  may  buckle  in  compression. 


2 . 2  PLANE  CRACK 

In  the  previous  section  2.1,  we  considered  the  effect  of  flaw 
volume  (porosity)  upon  bulk  modulus.  In  this  section  we  shall  show  how  to 
analyze  the  stress  field  around  a  plane  crack  in  an  infinite  plane  space 
subjected  to  simple  shear  at  infinity,  cf  sketch: 


within  the  framework  of  linear  theory,  the  problem  can  be  solved 
by  the  Kolosov-Muskhelishvili  technique  of  complex  variables.  It  Is  first 
convenient  to  map  the  Cartesian  frame  into  an  elliptic  frame  by  the  trans¬ 
formation: 


z  -  2  g  ) 


(42) 


where 

and 


C 

r.  - 

o  =  e 


(43) 

(44) 


Far  away  from  the  crack,  the  mapping  approaches: 

2  ^  p  e^^  (45) 

which  is  polar  in  nature.  For  this  problem  we  have  the  BC,  as  p  -* 


o  =  0  =  0,  T  =  S 
X  y  '  xy 


(46) 


These  stress  components  may  be  converted  by  tensor  transformation  into  the 
corresponding  polar  components ,  and  we  obtain 


as  P  ,  0  -  iT^  . -  iS  o 

P  P 


2 


(47) 


asP-O,  -  0 


(48) 


In  terms  of  the  two  complex  potentials  0  and  ip,  the  cartesian  displacements 
and  elliptic  stresses  are  given  by: 


2G(u  +  Iv)  =  (3-4v)<p-  -14,  -  0 

c  ^ 


o  +  o.=  2  (— 

P  '  2 


(49) 

(50) 


and 


0,  0; 

°p 


2  0--  ^rr 

-  ^  Iz  (— _ a 


C 


2 — ) 


(51) 
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I 


We  now  expand  <t>  and  ^  in  power  series  In  C  which  are  holomorphlc  for 
C  -♦  co;  we  then  equate  powers  of  a,  and  are  led  to  the  results; 


(52) 

(53) 


After  transforming  back  to  Cartesian  coordinates  we  obtain; 


xy 


(d+i)(c^-2^+i)  ^  f’  ({^+iKr:'*-2c^  +1) 

o  rk  ^  T  “ " "  ^ 


2  (£^  -  1)^ 


(C^ 


1)' 


2{V  -  1)^ 


r3  (C+^) 


C 


(C^  -  1)^ 


(54) 


What  we  desire  now  Is  an  average  value  of  t  within  some  unit  cell  of 

Xy 

outer  radius  R,  Inner  radius  1.  If  one  defines 


where 


<T  > 

S 

Jl  ^  V 

ir(R^  -  1). 


(55) 

(56) 


one  arrives  at  the  Incorrect  result; 


S 


(57) 


The  reason  that  <  turns  out  to  be  Independent  of  R  under  this  definition 
of  the  average  is  that  there  Is  a  stress  singularity  at  [p=  1;  ^  =  0,  ir].  This 
stress  singularity  (or  stress  concentration)  affects  the  reduced  stress  at 


&  =  j  ,  In  such  a  way  that  <t„„>  Is  a  constant.  In  order  to  define  <T  > 

^  xy  xy 

correctly,  one  has  to  treat  the  body  with  a  crack  In  it  as  a  "black  box"  and 
use  the  phenomenological  definition  of  shear  modulus,  namely; 


(58) 
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with 


<G  >  = 


(59) 


h  ^  ^eff  P 

ir  (R^  -  1) 


This  Integral  can  In  principle  be  done  but  the  work  Is  prohibitive .  In  the 
Interest  of  ‘'xpedlency,  we  turned  our  attention  to  another  problem  which 
can  be  done  in  a  more  straight  forward  fashion  with  equally  useful  results. 


2.3  EFFECTIVE  POISSON'S  RATIO  AND  YOUNG'S  MODULUS 

Consider  now  an  elastic  space  perforated  by  a  sphere  of  radius 
&  subjected  to  simple  tension  at  z  =  ±  in  the  absence  of  the  hole,  the 
solution  In  Cartesian  coordinate  Is,  of  course: 


In  spherical  coordinates,  the  solution  (r 

=  S  cs^  e 

—  S  sn  Q 

=  S  sn0  CS0 
°e»  0 

=  -  S  sn0  CS0 
u  =  p  r  [cs^0  -  vsn^0] 
V  =  -  sn0  CS0 
w  =  0 


(60) 

(61) 

(62) 

(63) 

w)  takes  the  form, 

(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 

(71) 
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0 


(72) 


At  the  hole ,  we  have: 

We  present  here  only  the  pertinent  details;  the  solution  Is  straight  forward. 
We  expand  u  and  v  In  spherical  harmonies,  and  write; 

00 

u  =  S  a  P  (cse)  (73) 

rt=0  "  " 

00 

V  =  L  PJ^  (csi»  (74) 

where  the  prime  denotes  differentiation  re  r_,  and  a and  depend  on  r_. 

We  use  the  method  of  equating  coefficients  of  P^^,  and  after  much  calculation, 
obtain: 


n  ,  4  -  5v 
^  2(7-5v) 


9 

2(7-5v) 


2 

X  cs  6  r  33  +  15v  ,a^\3 
'r  ^  7-  5v  ^  2  'r  ^ 


15  cs^e 
2  (7-5v) 


[7  (1)5  -  5  (1)3] 


(75) 


(4-15v)(H-  V) 
2  (7-5v) 


X  15  (1  +  V )  cs^Q  r  ^ 

2  (7-5 v  )  ^  'r  ^  ■  ^r  ^  J 

where  w  Is  the  Cartesian  displacement  component. 

The  effective  Young's  modulus  Is  defined  by: 

^eff  ^  ^  A+  B  cs^e  +  Ccs^e 

^  wE  (D  +  B^cs^e) 

where  {A,  B,  C,  D,  E^3  are  defined  above  In  (75)  and  (7  6). 
The  average  tensile  modulus  is  then  given  by; 


(76) 


(77) 
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<E> 

£ 


^  4«  BE  -  CD 

^  f 2  A+Bcs  9^Ccs  9  |  =  {^  +  - 

o  (D  +  E^cs^e)  r=b  *^^0  E^ 


+  [  - ^  ^  ^r-b 

/DE7  ^o  ^o  E“  ^o 
o  o 

2 

fA  f,  _ £  I  _o  ^  X  —  /JL  _ 2.  \  X  ^  1  /7Q^ 

-  ^D  "  3D  D  '3  ■  5D^  SD^r^b 

<  E  >  a  ^ 

Figure  1  shows  how  —  behaves  as  a  function  of  (^)  .  In 

<E> 

plotting  ,  it  was  sufficient  to  use  the  approximate  representation  given 
by  (79),  which  is  valid  for  0  ^  ^  the  limit  of  small  (^)^i  (81) 

behaves  like: 

“  -  1  2  i7  -,5v)  'b>  ‘®“) 

which  means  that  the  Initial  downward  slope  of  the  two  lines  shown  In 

Figure  1  for  (v  =  j  /  j)  and  respectively  \  ^  J  •  Equation  (80)  can  now  be 

used  In  a  fracture  analy  sis  as  long  as  .(|-)^  Is  small  with  respect  to  unity 

o  <E> 

In  a  sense  that  can  be  defined.  We  prefer  however,  to  represent  — p- 

^  3“ 

empirically  by  an  expression  that  fits  over  the  whole  range,  0  ^  (^)  ^1. 
Such  an  expression  Is  of  the  form: 

<|>  ~  1  .  .13  -  2y^l_5^  ,|,3  pj) 


27  -  12v-  15v 


Equation  (81)  assumes  that  the  two  curves  In  Figure  1  are  roughly  linear. 

The  expression  for  Poisson's  ratio  Is  obtained  In  the  following 
way.  The  Cartesian  displacement  is  given  by: 


U  V 


u  +  V  ctn  0 


u  -  vtn  0 
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80  that 


-V  =  ~  a  ^  VCtn  a 
eff”  ^  u  -  vtan© 


where  u  and  v  are  given  by: 


—  r  JL  .  1  “  2v  r  ,  o  _  < 

as  -  I  12  ^r^  +  -5+1^  a  ]+  P2  (cse)  (f-) 


5  (5  -  4v)  .a 


+  --Y'  (3.\  i  i  r  1 

^  6  (7  -  5v)  ^  3  a  ^ 

Substitution  Into  (84)  yields  after  averaging  at  r  =  b: 


where 


TTn  c  T  +  (— )  +  -  /a\^  .  1 

2  ^-5v)  ^  M7-5v)  ^b^  ■*’2  (88) 

«=  (|)^  (89) 

^  "  4  {7-5v)  t  ( b^  “  (b^  ]  (90^ 

Figure  2  shows  how  (87)  behaves  as  a  function  of  (-)^  In  nlr^tn 
(87)  It  was  sufficient  to  use  the  approximation: 


<V>  =  ^(l-i-2+lcf» 

B  '  3  B  ^  5  _2  '  "  ^ 

o 


In  the  limit  of  small  (|)  ,  (gi;  behaves  like: 


<V>  ^  V  -  v)  (14v  -  1 


2  (7-5v)  (p 


14 


which  means  that  the  initial  downward  slope  of  the  lines  in  Figure  2 
11  3  1 

(for  V  =  j)  are  respectively  ^  •  Over  the  whole  range  (91)  can  be 
better  represented  by  an  empirical  function  of  the  form: 


vci  V  -  k  (|) 


+  1 


ih 


(93) 


where  k  is  the  linear  coefficient  in  (92),  and  1  is  adjusted  so  that,  at 

(l)^=  1. 


<v>  = 


1  -f-  4v  -  5v  ^ 
9  -  4v-  Sv^ 


(94) 


which  is  the  value  obtained  from  (94). 


The  results  of  the  preceeding  analysis  can  be  summarized  by 
saying  that,  for  a  porous  material,  an  effective  Young's  modulus  and  an 
effective  Poisson's  ratio  can  be  defined  by  an  expression  of  the  form: 

2 

and  <v>  =  V  -  kj  Vp  +  k2  Vp  (96) 

where,  if  the  holes  Interact,  Vp  is  given  by: 

''p'=  <¥■>"  (87) 

with  b  and  ri  empirical.  Also  the  parameters  {k^  kj  k2  }  depend  on  v  as 
given  above . 


Before  using  <E>  and  <v>  to  determine  the  fracture  behavior  of 
a  randomly  flawed  material,  it  is  prudent  to  point  out  that  the  representa¬ 
tions  (95)  and  (96)  are,  strictly  speaking,  only  valid  in  a  simple  tensile 
field.  It  would  be  necessary  to  determine  these  same  parameters  for  two 
or  three  other  types  of  loading  and  to  strike  some  sort  of  an  average  behavior. 
In  general,  however  the  algebra  Involved  is  cumersome  and  tedious  and,  at 
the  best,  we  must  think  of  k^,  kj^,  and  k2  as  seml-emplrlcal  parameters 


anyway.  Thus  we  prefer  to  eschew  the  analysis  of  other  model  flaw 
problems  and  will  go  ahead  on  the  basis  that,  roughly  speaking,  effective 
Young's  Modulus  decreases  linearly  with  flaw  content,  and  effective 
Poisson's  ratio  decreases  parabollcally  with  flaw  content. 

2.  4  INCORPORATION  OF  EFFECTIVE  MODULI  IN  THE  GRIFFITH 

.  CRITERION 

Again,  we  start  slowly  and  look  at  a  model  problem,  one  which 
was  first  treated  by  Max  Williams^,  and  v^'hlch  is  the  first  problem  con¬ 
sidered  in  this  section,  namely,  the  hydrostatic  tension  of  a  sphere. 

Imagine  now  that  as  the  sphere  is  stretched  radially  by  a  monotonlcally 
increasing  stress  at  its  outer  boundary  (r  =  b),  that  the  inner  boundary  (r=a) 
starts  to  ablate  (fracture)  when  some  critical  stress  concentration  is  reached. 
In  this  analysis,  we  ignore  inertia  and  do  a  quasi-static  analysis  (the 
algebraic  details  are  suppressed). 

In  the  case  of  linear  theory,  it  is  easy  to  show  that,  as  the  stress 
is  increased  monotonlcally  from  zero,  a  critical  stress  is  reached  which 
depends  on  the  initial  flaw  size,  and  which  is  determined  by  the  Griffith 
criterion  to  have  the  form: 

(98) 

where  a^  is  the  critical  hole  size  and  F is  the  specific  surface  energy. 

At  this  point  the  linear  theory  predicts  that,  in  order  for  a^  to 
increase,  S  must  be  decreased;  and  conversely,  as  S  is  increased,  a^ 
decreases.  Thus  (98  )  possibly  holds  only  at  the  point  of  initiation  of 
flaw  growth,  and  predicts,  as  is  known  from  studies  on  glass,  that  the 
critical  stress  depends  inversely  on  the  square  root  of  the  hole  size.  In 
order  to  find  how  the  hole  grows  as  S  increases  indefinitely,  one  must  do 
a  geometrically  exact  analysis.  In  order  to  gain  insight  into  how  the  hole 
grows,  we  present  here  an  exact  geometrically  non-linear  analysis  for  an 
incompressible  neo-Hookean  material.  Following  the  same  approach  as 
used  above  in  Section  2.1,  it  is  easy  to  show  that  the  Griffith  Criterion 
predicts  that: 
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(99) 


1/3  r  ,  1  2/3  \l-l  -4/3 

V  +  (1+-V-)  -3] 

p  ''p 


where 


P  b  '^p 


Xl-  l+3c^ 


K  -  1  -1/3 

V-)  ] 


(100) 

(101) 


Ofi 

In  the  limit  of  small  strain,  l.e. ,  as  — -  -*  0  we  have: 

''p 


r  5/6 

S=  ^Eb-  ^p 


(102) 


s  =  bJIZ 

G  “  ^  Eb  “l/e 


(103) 


which  relations  are  valid  for: 


V  >> 

'Gb' 


c. 

In  the  limit  of  larQe  strain,  l.e.,  as  3  —  co  \y@  have: 


(104) 


1  2r  ^^2  _ 

*b"  3  ^ 


(105) 


I  - 


(106) 


Which  relations  are  valid  for: 

2r  3 


(107) 


Thus,  in  summary,  for  large  strain,  Vp  increases  quadratlcally  v/lth  strain, 
and  for  small  strain,  the  rate  of  increase  drops  to  the  (6/5)^^  power.  This 
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behavior  Is  sho^/vn  In  Figures  3  and  4  where  v  Is  plotted  versus  log  for 
various  values  of  (4r/Gb).  In  Figure  5  we  present  a  crossplot  of  log  c 
(Vp-  1)  versus  by  (4r/Gb).  In  the  smaU  strain  Hml  ,  (103)  compares  favor¬ 
ably  with  (  98)  within  a  factor  of  two,  whereas  behavior  beyond  20%  strain 
Is  Just  not  forthcoming  from  the  linear  theory.  Whether  or  not  one  uses  the 
nonlinear  theory  depends  on  the  value  of  (F/Gb), 

In  general,  for  weakly  cohesive  materials,  linear  theory  Is 

adequate. 

Let  us  now  redo  this  problem  in  the  framework  of  linear  theory, 

and  replace  the  cavltated  ^hell  by  a  solid  shell  of  radius  b  with  an  effective 
bulk  modulus  given  by  (1). 

A  simple  calculation  shows  that; 


(108) 


Thus  a  Griffith  calculation  based  on  effective  elastic  constants  rather  than 

flaw  geometry  leads  to  a  solution  of  the  same  form  as  that  obtained  In  the 

"flaw”  calculation,  except  that  the  length  scale  a^  Is  replaced  by  b.  In 

effect,  this  means  It  Is  equivalent  to  filling  up  the  void  with  material  of  the 

ooireot  effective  modulus.  Thus  the  randomly  flawed  theory  leads  to  reason 
able  results. 

It  Is  our  Intention  to  examine  In  similar  fashion  the  problem  of 
simple  tension  of  the  space  with  a  hole  replaced  by  "effective"  material. 

If  this  checks  out,  we  will  then  program  the  relations  (  95)  and  (  9)6)  Into 
the  SHEP  code,  and  will  then  run  several  wave  problems  using  the  Griffith 
criterion  for  materials  with  simple  geometry  and  moduli  which  depend  on  flaw 
fracture.  It  Is  then  our  Intention  to  define  a  critical  a  or  v  at  which 
sufficient  damage  is  assumed  to  have  occurred  so  that  the  fractured  material 
IG  now  ready  for  excavation  by  shoveHng  processes.  Thus  we  have  here  a 
simple  damage  theory,  a  precursor  to  excavation. 
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The  present  theory  has  so  far  been  developed  for  an  Isotropic 
elastic  continuum.  In  the  forthcoming  period ^  we  propose  to  examine  both 
an  orthotropic  elastic  material  as  well  as  an  elasto-plastlc  material. 


3.  TECmnCAL  PROGRESS:  DISCRETELY  FAULTED  MKDIA 


In  the  work  completed  on  Contract  No.  H0210019,  It  was  shown 
that  a  computer  solution  obtained  by  SHEP  provided  excellent  agreement  with 
the  same  features  of  the  analytical  solution  derived  by  Achenbach  and  later 
extended  by  Blatz.^  One  feature  of  the  computer  solution  Is  not,  however, 
desirable.  Since  the  SHEP  code  Is  a  difference  equation  representation  of 
the  continuum  equation.  It  Is  naturally  resonant  and  “coughs  back"  both 
oscillations  In  the  wave  profile  as  well  as  spreading  of  the  wave  profile. 

This  behavior  Is  shown  In  Figure  6  where  the  computer  solution  Is  compared 
with  the  theoretical  solution  for  the  velocity  profile  from  the  crack  tip  to 
the  bow  wave  front.  Such  a  numerical  behavior  makes  It  Impossible  for  us 
to  compute  precisely  enough  the  Integrals  In  the  algorithm  for  crack  propa¬ 
gation.  An  attempt  to  smooth  out  the  curve  is  being  made  with  the  aid  of 
cubic  viscosity,  an  extension  of  the  form  Introduced  by  Von  Neumann  and 
Rlchtmeyer. 

If  we  decompose  the  stress  tensor  Into  elastic  and  viscous  com-  ■ 
ponents,  we  find: 


(109) 


where  t  Is  the  dissipative  stress.  This  component  generally  can  be 
written  as; 

^  ^  ^  _  (110) 

where  d  is  the  strain  rate  tensor,  and; 


«k=  ftk  ^id' “d^  iiy  ^  2  (111) 

with  {l^,  II^,  IIIjj}  the  three  Invarients  of  the  strain  rate  tensor.  Von  Neu¬ 
mann  and  Rlchtmeyer  took; 
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s 


Figure  6.  Computer  Oscillations  for  Out  of  Plane  Motion. 
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(112) 


*D”  2 '’o'*  (i-iy 


where  is  a  kinematic  viscosity  and  T  Is  a  volumetric  relaxation  time. 


We  take  tj,  =  2v^  p  (1-fiU^)  d 


(113) 


where  fi  Is  a  higher  order  rate  constant. 


Eguatlon  (112)  Involves  quadratic  and  linear  terms  /  (113)  Involves  cubic  and 
linear  terms.  Equantlon  (113)  provlces  a  twofold  advantage  over  (112).  First 
It  can  be  used  for  Incompressible  or  Isochorlc  deformation  (such  as  our  model 
antiplane  shear  problem).  Equation  (112)  can  only  be  used  for  compressible 
media.  Equation  (113)  contributes  to  both  the  Inplane  and  out-of-plane 
terms.  In  (112)  the  quadratic  term  contributes  only  to  the  Inplane  terms. 

Thus  (113)  fills  the  bill  for  our  problem. 


The  coefficients  and/3  are  given  by: 


At 


2 


(114) 

(115) 


Equations  (113)#  (114)#  and  (115)  have  been  programmed  and  an  optimization 
study  of  the  constants  and  C2  Is  underway. 
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